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Abstract 

We have shown that non-zero U e s is generated in a see-saw type softly broken A4 symmetric 
model through a single parameter perturbation in mjj in a single clement. We have explored 
all possible 9 cases to study the neutrino mixing angles considering the best fitted values of 
Atoq and AmJ tm with all parameters real. We have extended our analysis for the complex 
case and demonstrated large low energy CP violation (Jcp — 10~ 2 ) and m ee in addition to 
mixing and mass pattern. We have also investigated leptogenesis and for a reasonable choice 
of model parameters compatible with low energy data, WMAP value of baryon asymmetry 
6 x 1CP 10 is obtained for right handed neutrino mass scale M ~ 10 13 GeV. We have obtained 
a relation among the phases responsible for leptogenesis and have shown that those phases also 
have correlations with low energy CP violating phases. 

PACS number(s): 14.60.Pq, 11.30.Hv, 98.80.Cq 

1 Introduction 

In recent time people are too much interested to find some flavor symmetry in order to generate 
mass and mixing pattern of fermions. Continuous symmetry like L e — — L T [1 , — L T [2] 
symmetry and most popular discrete symmetry, fi — r exchange symmetry (S^ 1 ")!^ have got some 
success to describe mass and mixing pattern in leptonic sector. To avoid mass degeneracy of \i and r 
under symmetry, E. Ma and G. Rajasekaran in [5] have introduced first time the A<± symmetry. 
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After this paper, a lot of work have done with this symmetry j5]-jll]. After introduction of tri-bi 
maximal mixing pattern (sin #12 = l/\/3, sin #23 = — 1/ \/2, sin #13 = 0)0, people have tried to 
fit this mixing pattern through the A4 symmetry. In an well motivated extension of the Standard 
model through the inclusion of A4 discrete symmetry tri-bi maximal mixing pattern comes out 
in a natural way in the work of Altarelli and Feruglio [6j. More precisely, the leptonic mixing 
arises solely from the neutrino sector since the charged lepton mass matrix is diagonal. The model 
[6J also admits hierarchical masses of the three charged leptons whereas the neutrino masses are 
quasi-degenerate or hierarchical. Although the model gives rise to #13 = (f/ e 3 = 0) which is 
consistent with the CHOOZ-Palo Verde experimental upper bound (#13 < 12° at 3a), however, the 
non-zero and complex value of U e 3 leads to the possibility to explore CP violation in the leptonic 
sector which is the main goal of many future short and long baseline experiments. Within the 
framework of SU{2)l x U(1)y x -A4 model, non-zero U e 3 is generated either through the radiative 
correction [21] or due to the introduction of higher dimensional mass terms [6] . Generation of non 
zero complex U e s and possibility of non-zero CP violation has been extensively studied in [T7] for 
the proposed model of Altarelli- Feruglio |6j with explicit soft breaking of A4 symmetry |21j . 

In the model [6J the authors showed that the tri-bi maximal mixing pattern is also generated 
naturally in the framework of see-saw mechanism with SU{2)l x U{1)y x symmetry. Exact 
tri-bi maximal pattern forbids at low energy CP violation in leptonic sector. The textures of mass 
matrices in [6J could not generate lepton asymmetry also. In the present work, we investigate the 
generation of non-zero C/ e 3 through see saw mechanism by considering a small perturbation in mo , 
the Dirac neutrino mass matrix, keeping the same texture of the right-handed Majorana neutrino 
mass matrix as proposed in Ref.|6 . At first, we have studied in detail perturbation of mo by adding 
a small parameter at different entries of mo and see the variations of three mixing angles in terms 
of other model parameters considering all of them real. We extend our analysis to the complex case 
for a suitable texture. We study detailed phenomenology of neutrino mass and mixing including 
CP violation at low energy, neutrinoless double beta decay and leptogenesis. Our approach to get 
nonzero U e % is minimal as we break A4 symmetry explicitly by single parameter in single element 
of me. Generation of CP violation at low energy as well as high energy is also minimal as we 
consider only one parameter complex. 



2 The Model of Altarelli Feruglio with See-Saw: light neutrino 
phenomenology and leptogenesis 

We consider the model proposed in jB], which gives rise to diagonal and M\ (the charged 
lepton mass matrix) along with a competent texture of Mr and after see-saw mechanism and 
diagonalisation gives rise to tri-bimaximal mixing pattern. The model consists of several scalar 
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Table 1: List of fermion and scalar fields used in this model, I = e, fi, t. 



fields to generate required vacuum alignment to obtain tri-bimaximal mixing. In Table I., we have 
listed the scalar fields and their VEV's and representation content under all those symmetries. 
The model is fabricated in such a way that after spontaneous breaking of A4 symmetry, the S% r 
symmetry remains on the neutrino sector and the charged lepton sector is invariant under Z3 
symmetry. Consider the Lagrangian of the model |6j, 

£ = ^(<hfi L )e R h d + y ^{<fri> l L )' m h d + ^^t^lY'trK + ffi L N R h u 

+x A ^N c L N R + x B( PsN c L N R + h.c. (2.1) 

After spontaneous symmetry breaking, the charged lepton mass matrix comes out diagonal with 
m e = Ve ^^ d 1 "fy = Vti J^^ d 1 an d itt-t = Vr ^pi^ ■ The neutrino sector gives rise to the following Dirac 
and Majorana matrices 

-D/3 -D/3 

m D = f^L\Q 10] M R =\ -D/3 2D/3 A-D/3\ (2.2) 





A — D/3 2D/3 

where A = 2xau, D = 2x R vs- The structure of light neutrino mass matrix can be obtained from 
see-saw formula: 

/ 2(D+A) u u 

\l~h (2.3) 





M v = -moM^m^ = U TB 



2{D-A) 
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where, 



Utb 



V 



(2.4) 



This is clear from Eq ]2.3| that Utb is the diagonalising matrix for light neutrino mass matrix M v . 
The form of Utb is i n Eq ]2.4| which is nothing but the so called tribimaximal mixing matrix. From 
]2.3| we have the eigenvalues of M u : 



mi 



2(D + A) 



m 2 



2A 



m 3 



2(D — A) 



(2.5) 



From Eq. 2.4 we have the mixing angles sin #12 = l/v3, sin 623 = — l/v2 and sin #13 = and from 



2.5| we get the solar and atmospheric mass squared differences as 

m 2 k{k + 2) 



Am: 







mj 



ml 



Am 



atm 



ml 



(1 + k) 2 
m\k(2 - k) 
(1-k) 2 



(2.6) 



p2 2 

where D = kA, mo = an d all parameters are real. From the experiments we know Amg is 
positive and dictates either k > or k < —2. If k > 0, then it has to be small in order to generate 
small value of Am|, provided m\ is not too small as Anig. But small positive k corresponds to 
same order of magnitude of Am|, and Am\ tm which is not acceptable according to the experimental 



results. Now k > only acceptable for m. 



2 atm obtained 



— Attiq and hierarchy of Attiq and Am, 

near k ~ 1. This corresponds to normal hierarchical 



2.6 



with the singular nature of Am^ im as in Eq. 
mass spectrum. Again for m\ 3> Am|,, k < —2 is the physical region. This region of k makes 
Am\ tm < which is so called inverted ordering of neutrino mass pattern. Again k + 2 should take 



small value in order to generate small value of Amg. For one complex parameter D 
can write the mass differences in the following form 



De l 



wc 



9 mikik + 2 cos 6 

Am„ = u , \ : - 

l + /c 2 + 2fecoS( 



Am 



atm 



m.0&;(2 cos <p — k) 
1 + k 2 — 2k cos 6 



(2.7) 



In the complex case, positivity of Awig can be obtained either with k > and cos<f> > —k/2 or 
with k < and cos (f> < —k/2. For the first case with m,g ~ Am|, and with cos <j) ~ (1 + k 2 ) /2k one 
can have normal hierarchical mass spectrum. For the second case hierarchy will be inverted and 
k + 2cos(t!> have to be small. In both case k should take the value such that the 1 > cos 6 > — 1 



range also satisfy. The mixing pattern is tri-bi maximal Eq,2.4 and it is independent to the fact 
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whether the parameters are real or complex. In this mixing pattern U e s = and non-zero complex 
U e 3 is a basic requirement to see the non-zero Dirac CP violation. 

Now we concentrate on the issue of leptogenesis in this model. The decay of right handed heavy 
Majorana neutrinos to lepton(charged or neutral) and scalar (charged or neutral) generate non-zero 
lepton asymmetry if i) C and CP are violated, ii)lepton number is violated and iii) decay of right 
handed neutrinos are out of equilibrium. We are in the energy scale where A4 symmetry is broken 
but the SM gauge group remains unbroken. So, the higgs scalars both charged and neutral are 
physical. The CP asymmetry of decay is characterized by a parameter £j which is defined as 



Ni->l+<f>- ,if<jfl* 



1 N^l-<j>+,v,<t> + ^Ni-+l+^>-,vf<jfl* 

Spontaneous A4 symmetry breaking generates right handed neutrino mass and the mass matrix 



(2i 



Mr obtained is shown in Eq. 2.2 We need to diagonalize Mr in order to go into the physical basis 
(mass basis) of right handed neutrino. This form of Mr gives the diagonalising matrix in the tri-bi 



maximal form Utb m Eq. 2.4 



U TB M R U TB = diag(Mi, M 2 , M3) = diag(^4 + D, A, D-A), 



(2.9) 



however, the eigenvalues are not real. We need to multiply one diagonal phase matrix Up with 
Utb- Hence, diagonalising matrix V = UtbUp relates the flavor basis to eigen basis of right handed 
neutrino: 



N lR = J2V l *N m . 



i=i 



In this basis the couplings of Nr with leptons and scalars are modified and it will be: 

m' d = m d V*. 



(2.10) 



2.11 



At the tree level there there are no asymmetry in the decay of right handed neutrinos. Due to the 
interference between tree level and one loop level diagrams, the asymmetry is generated. There 
are vertex diagram and self energy diagram to contribute to the asymmetry |47[ 48 . The vertex 
contribution is : 



-V 



1 



ij > 



fxij { 1 - (1 +Xy)ln(l + 



1 



(2.12) 



and the self energy part is : 



ef = 4W^£ /m(/ ^ )X 1 



(2.13) 
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where x it = M 2 /M 2 and 



h 



(2.14) 



The key matrix, whose elements are necessary to calculate leptogenesis, is h. In this model mr> is 
diagonal and proportional to identity. Hence, h matrix is real diagonal and it is also proportional to 
identity matrix and it is independent of the form of V. The terms for decay asymmetry generated 
by "i" th generation of right handed neutrino Nm for both vertex and self energy contributions 

and Eq. 



are proportional to Im(/i 2 ) (where j ^ i) as in Eq 



2.12 



2.13 



All off-diagonal elements 



of h are zero. So, decay of all three generation of right handed Majorana neutrinos could not 
generate lepton asymmetry. So, in this model of A\ symmetry tri-bi maximal mixing pattern is not 
compatible with the low energy Dirac CP violation as well as high energy CP violation. In order 
to obtain non-zero #13, low energy Dirac CP violation and leptogenesis we need to break the A4 
symmetry through not only spontaneously but also explicitly introducing some soft A± symmetry 
breaking (soft in the sense that the breaking parameter is small to consider A4 as an approximate 
symmetry) terms in the Lagrangian. 



3 Explicit At symmetry breaking and real parameter analysis 

We consider minimal breaking of A4 symmetry through a single parameter in a single element of 
mo keeping Mr unaltered as 

^A 4 breaking = f ^qL^^R^u- (3.1) 

We introduce the breaking by small dimensionless parameter e to the (a,/?) element of Dirac type 
Yukawa term for neutrino. After spontaneous SU(2)l x U(1)y symmetry breaking it modifies only 
one element (a,/3) of tbd of neutrino. There are nine possibilities to incorporate the breaking 
parameter e in mo- We know that after spontaneous A4 symmetry breaking, a residual SrT 
symmetry appears in neutrino sector. There is a special feature of SrT symmetry which ensures 
one and one maximal 7r/4 mixing angles. There is one task to check whether our newly introduced 
explicit breaking term can break S% T symmetry or not. This is important because we need non-zero 
#13. We have seen that in one case out of the nine possibilities, residual symmetry remains 
invariant. This is af3 = 11 case. In other cases S% T symmetry is broken and one expect non-zero #13 
from those cases. Primarily, we consider that all parameters are real. We want to study the mixing 
pattern and want to see its deviation from tri-bi maximal pattern considering experimental value of 
mass squared differences of neutrinos. We explore all nine cases including afi = 11 case. Although 
a(3 = 11 case could not generate non-zero #13, however, we want to see whether this breaking can 
reduce the tri-bi maximal value of #12 (— 35.26°) to its best fit value (~ 34°) or not along with the 
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special feature 9 13 = and #23 = tt/4. Here, we explicitly demonstrate the procedure for a single 
case and for the other cases expressions for eigenvalues and mixing angles are given in Apendix. 

(i) Breaking at '22' element : In this case, the structure of mo is given by 



m D 




(3.2) 



and after implementation of see-saw mechanism keeping the same texture of Mr, three light neu- 
trino mass eigenvalues come out as 



mi 



m 2 



2(D + A) 

and the three mixing angles come out as 
1 e(2A + D) 



1 + 



2c 



m 3 



sin 0i2 = 



V3 



3V3D 
sin 6>i3 



sin 6*23 = - 
D 



if 

3\^2A 4A-2D 



1 

V2D 



D 



2(D — A) 

V2D 



(1 + 6) (3.3) 



+ 



3 1 2V2A 4A - 2D 



(3.4) 



Assuming a relationship between the parameters D and A as D = kA we rewrite in a convenient 
way the above three mixing angles as 



sin 0i2 = 



1 



e(2 + k) 



3VSk 

and the mass-squared differences are 



sin fc>23 = - 



1 | ek(4 - k) 



V2 6^/2(2 - k) 



sin 9 



13 



ek{l - k) 
3\/2(2 - k) 



(3.5) 



Am; 



m. 











3(1 + A;) 2 



Ami 



m 







3k(k + 2) + 2e(2k 2 + 4k + 1) 
3A:(2 - fc) - 2e(2A; 2 - 4k - 1) 



w 3 (! _ fc )2 

where mo = / 2 f 2 /2^4. Defining the ratio R in terms of mass-squared differences we get 
R 



Am| 



_ (fc - l) 2 [3k(k + 2) + 2e(2fc 2 + 4k + 1)] 

Am 2 tm " (k + l) 2 [3A;(2 - k) - 2e(2k 2 -4k- 1)] 

which in turn determines the parameter e as 

3fc [«(2 -£:)(£; + l) 2 - (A; -l) 2 (fc + 2)] 



e = 



(3.6) 



(3.7) 



(3.8) 



2 p - l) 2 (2fc 2 + 4k + 1) + + 1) 2 (2A; 2 -4k- 1)] 

Similarly, we have evaluated all other possible cases which we have listed in the Appendix. 

Now the e parameter is determined in terms of R and k and we substitute it to the expressions 
for mixing angles. Thus it is possible to explore all three mixing angles $12, 023 and $13 in terms 
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Figure 1: Plot of 9 12 with respect to k . We keep Am| 2 and AtttJi to their best fit values. 
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Figure 2: Plot of #23 with respect to k. We keep Am| and Am^ to their best fit values. 
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Figure 3: Plot of #13 with respect to k. We keep Am| 2 and Am^ to their best fit values. 

of R and k. Particularly, the deviation from tri-bimaximal mixing depends only on R and k. For 
the best fit values of the solar and atmospheric mass squared differences (R ~ 4 x 10 -2 ), we have 
shown in Fig. [T] Fig. [2] and Fig. [3] the variations of #12, #23 , #13 verses k, respectively. We have 
studied all nine possible cases and shown in the plots. First of all, non-zero value of #13 is obtained 
if we allow A4 symmetry breaking terms explicitly in any one of the '12, 13, 21, 31, 22, 33' element 
of the Dirac neutrino mass matrix and those are the cases of our interest. In the present analysis, 
we have shown that non-zero #13 is generated in a softly broken symmetric model which leads 
to deviation from the tri-bimaximal mixing. In A4 symmetric model #13 is zero because of residual 
symmetry in neutrino sector after spontaneous breaking of A^ symmetry. Apart from the 
explicit breaking of A4 symmetry at 11 element, is broken for all '12, 13, 21, 31, 22, 33, 23, 
32' cases. Furthermore, perturbation around '23', '32' elements also lead to zero value of #13 at the 
leading order although Si^ symmetry is broken, non-zero value is generated if we consider higher 
order terms of e 2 which are too tiny and hence, discarded from our analysis. We include 11 case 
for completeness which preserves symmetry and hence generates #13 = and #23 = tt/4. It 
only shifts #12 from the tri bimaximal value, but it cannot be able to go towards the best fit value 
of solar angle , 34° . 

If k = —2, then we get e oc R, and thereby, the value of #13 is very small also #12 will hit 
the exact tri-bimaximal value in some cases. The effect of variation on the mixing angles around 
k = —2 are asymmetric. For some cases (for example 23, 32) #23 changes very fast in the k > —2 
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region. So, we explore the mixing angles with the range —2.15 < k < —1.95. We choose the most 
feasible cases in which perturbation is applied around '12', '13' elements, because in those cases, 
variation of k encompasses the best-fit values of #12 and 623- Although, in the '21', '31' cases, the 
value of 623 touches the best-fit value 7r/4, however, 612 far apart from the best-fit value. In order 
to achieve large #13, we have to choose the '21', '31' cases, but we have to allow the variation of 9\2 
around as large as 37.2°. In case of '22', '33', the structure of m£> is still diagonal and also we can 
get larger #i3(upto 6°) and also 612 is within ler(36°), however, 623 will reach 3cr(54°) value. 

In summary, we have shown that non-zero #13 is generated in a symmetric model which leads 
to deviation from the 'tri-bimaximal' mixing through see-saw mechanism due to the incorporation 
of an explicit symmetry breaking term in mj). The breaking is incorporated through a single 
parameter e and we have investigated the effect of such breaking term in all nine elements of ma. 
Some of them generates still zero value of #13 and rest of the others generated non-zero #13. We 
expressed all three mixing angles in terms of one model parameter and showed the variation of all 
three mixing angles with the model parameter k. We find breaking through '12' and '13' elements 
of m/) are most feasible in view of recent neutrino experimental results. 



4 Complex extension: Light neutrino phenomenology and Lepto- 
genesis 

In this section, we consider one of the parameter is complex and out of all nine cases as mentioned 
earlier, we investigate one suitable case arises due to '13' element perturbation. This is one of 
the suitable positions of breaking justified from real analysis. Again this extension is minimal to 
generate non-zero CP violation because we consider only one parameter complex. We take D as 
complex: D = De 1 ^. Hence, the form of mn and Mr under explicit A4 symmetry breaking with 
complex extension are : 



m D = f 



/I 





e\ 







1 





Mr 


\o 





J 





(A+ (2De < *)/3 
-(De**)/3 
-(De**)/3 



-{Dj*)/Z -(De^)/3 \ 
(2De^)/3 A-{De^)/Z 
A-(De^)/3 (2De^)/3 ) 



(4.1) 



4.1 Light neutrino phenomenology 

Using the see-saw mechanism we get the light neutrino mass matrix as 



R 
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Utb 



V 



1 _ 2e 2De i <t'(2A+De i 't')e 2 

De't+A 3(£»e**+A) 9A(A 2 -Z) 2 e i2 </>) 
(A+2De i <t>)e _ y / 2De^(2A+De^)£ 2 
3V2(De^+A) 9A(A 2 -D 2 e i2 <t>) 



{A+2De x <t>)e 
3v / 2(De^+A) 

1 _L 2£ 

A 3A 



_ v / 2De"»(2A+De^)£ 2 

9A(A 2 -£> 2 e i2 *) 
£>e^(2A+£>e^)<; 2 
9A(A 2 -D 2 e i2r l>) 
— e 





V6(A-De i <l>) 
1 



/7 T 



V3(A-De^) ^(A-De^) D&9—A ' 

We need to diagonalize the M v to obtain the masses and mixing angles. The eigenvalues are same 
as we have in the real case and only difference is that the D is complex now. We explicitly write 
down the complex phase in the mass matrix. The obtained eigenvalues are : 

fvj (. 2e\ _«/ 2e\ fvj 

~2{De^ + A)\ 3) m2 ~ 2A\ + 3) ™ 3 ~ 2{De^ - A) 

where we keep terms upto first order in e. Now with D = kA, mo = f 2 v 2 /2A and keeping term 
upto first order in e we get the three neutrino mass squared as 



mi 



(4.3) 



mi 



mg(9 - 12e) 



9(1 + k 2 + 2fccos <j>) 



JB2 



mg(9 + 12e) 



\m3\ 



m\ 



1 + k 2 - 2/ccos, 



(4.4) 



Using those expressions we get the mass squared differences and their ratio which are, 

m 2 {9k(k + 2 cos (j)) + 12e(2 + k 2 + 2k cos <p)} 

9(1 + k 2 + 2/ccos 4>) 
m§{9fc(2 cos (/) — k) - 12e(l + /c 2 - 2/c cos 0)} 



Am| 



™2 



Am 2 



aim 



rn 3 | - m 2 | 



and 



i? 



Am 2 

Awiatm 



1 + fc 2 - 2kcoscj) 
l + k 2 + 2k cos (j) 



9(1 + fc 2 - 2fccos0) 

9fc(fc + 2 cos 0) + 12e(2 + /c 2 + 2/c cos . 

12e(l + fc 2 - 2/ccos. 



(4.5) 



(4.6) 



9fc(2 cos </> — k) 

The mixing angles are obtained from diagonalisation of M v . We solve the equations of the form 
M u \ mi)* = mi | mi ). These | mj ) will give the columns of the diagonalising unitary matrix U. 
Throughout our calculation we assume that breaking parameter e is small. We have the nonzero 
U13 which is proportional to e. So, the values of Uu and U23 will give the solar and atmospheric 
mixing angles, respectively. The expressions for the mixing angles come out as 



sin U12 



\U- 



1 



12 



+ 



3^3 



2 + Ak 2 + 2k 4 + 9k cos <j) + 10fc 2 cos 2 (f> + 9k 3 cos (f> 
k(l + k 2 + 2k cos 4>){k + 2 cos </>) 



(4.7) 



sm #23 



23 



+ 



k 3 + k + 2k cos 2 — 3/c 2 cos (f> — cos 
\f2 ' 6\/2cos^ ' "fc 3 + k + 4/ccos 2 </> - 4/c 2 cos0 - 2cos<; 

13 ^ 13 ^ 3\/2 cos 0(/c 3 + k + 4fccos 2 — 4fc 2 cos <p — 2 cos <fi) 
(k 4 + k 2 - 3k 3 cos </> - k cos + 6k cos 3 0-3 cos 2 - k 2 cos 2 </>) 2 
2 



+ sm 



3 + k + 6k cos 2 (f) — 5k 2 cos <f> — 3 cos </>) 2 



(4.8) 



(4.9) 
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Figure 4: Plot of #12 with respect to k and eft for the breaking of in '13' element of mp. We 
keep Am§ 2 and Am^ to their best fit values. 



From the expression of mixing angles it is clear that the deviations from tri-bi maximal are first 
order in e. The independent parameters in this model are A, D, /, e and (ft. Alternatively the 
independent parameters are Mq, mo, k, e and (ft (where Mq = A,ijiq = f 2 v 2 /2A, k = D/A). In 
the above analysis of light neutrino mass and mixing, scale Mq did not appear explicitly. We have 
four well measured observable which are Am§,, Am„ (m , #12 and #23, and, thus, in principle it is 
possible to determine four parameters mo, k, e and (ft an d we are able to predict the other less 
known observable such as angle #13, CP violating parameter Jqp etc. It is difficult to get inverse 



relations of those observable. From the expression of R in Eq. 4.6 we easily obtain the expression 
for e as 



3/c [R(2 cose/) - fe)(l + k 2 + 2kcos(f)) - (2cos0 + k)(l + k 2 - 2fccos</>)] 
4 [R(l + k 2 - 2k cos <f>) (1 + k 2 + 2k cos </>) + (2 + k 2 + 2k cos 0)(1 + k 2 - 2k cos . 



(4.10) 



Now using the relation of Am 2 tm with the parameters Eq. 



4.5 



we get the expression for m\: 



m 2 



9(1 + k 2 - 2A;cos 4>)Am 



atm 



9k{2 cos (p-k) - 12e(l + k 2 - 2k cos . 



(4.11) 



where e is in the form of Eq. 



4.10 



Thus, e and rrin depend on the parameters k, cj) and experimentally 



known R. Extraction of k and (ft from other two known mixing angles is little bit difficult. Rather 
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Figure 5: Plot of #23 with respect to k and (f> for the breaking of in '13' element of mp. We 
keep Am§ 2 and Am|i to their best fit values. 




Figure 6: Plot of #13 with respect to k and <j> for the breaking of A4 in "13" element of mr>. We 
keep Am,32 and Am|i to their best fit values. 
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we have plotted #12 and #23 with respect to k and cf> and obtain the restriction on the parameter 



space of k and tfi. From the expression of sin #12 in Eq. 4.7 we are seeing that there is a factor 
k + 2cos0 in the denominator. For k > — 2 there will be a <j> for which the quantity k + 2cos0 
becomes zero. Hence, we should keep k < —2. Again the factor k + 2 cos (j) should be small to ensure 
that e is also small. It justifies our whole analysis because we consider first order perturbation as 
we considered symmetry A4 remains approximate. We consider the range —2.1 < k < —2.0 and 
< <p < 22°. From Fig. [4] we see that #12 changes from the tri-bi maximal value 35.6° to 31°. 
Near <f> = 10° it crosses the best fit value 34°. We have plotted #23 in Fig. [5] The variation of #23 
is from 47° to 44.8° for the same range of k and (j). The best fit value 45° is remain within range 
of variation and it is in the low k low <f> region. The plot of #13 is in Fig. |6] Value of #13 remains 
within 2° for the same range of k and (f> and the model predicts #13 is very small but non-zero. 
Question may arise whether such small value of #13 can generate observable CP violation or not. 

Keeping all those constraints in view next we explore the parameter space of CP violation 
parameter Jcp- The parameter Jcp defined as 



[ Im\h! h! h! 1 

J CP = - sin 29 12 sin 2023 sin 20 13 cos 6 13 sin 5 CP = . 2 L , 12 f. 3lJ 2 (4.12) 

where b! = M^M^, 5c p is Dirac phase. This Jqp is associated with CP violation in neutrino 
oscillation and is directly related to Dirac phase of mixing matrix. The from Eq. (4.3) we can 
express M v matrix in terms of e and mo and k and 0. With the expressions for e and mo we can 
easily obtain the M v and hence b! = M U M^ completely in terms of k and <f>. Hence, similar to the 
mixing angles Jcp will be also only function of k and (j>. We have plotted Jcp in Fig. [7] where 
the values are normalized by a factor 10 -3 . For the same range of k and cf> the model predicts 
Jcp up to the order of 10~ 2 which is appreciable to observe through the forthcoming experiments. 
Inverting the expression of Jcp, the phase 5c p is extracted in terms of k and cj) and it is plotted in 
Fig. [8j We see that the value of 5c p is large upto 90° and, thereby, compensates small 6*13 effect in 
Jcp and makes it observable size. One important discussion we have to make about the range of 
(j) and k. One can ask why we are keeping ourselves small range of those parameters where larger 
(j) can enhance the #13 and Jcp- We have studied that the larger value of <p and also k in negative 
become responsible for breaking the analytic bound sin 5c p < 1- So, we keep ourselves in shrinked 
parameter space which keep 612 and #23 in exceptionally good values according to the experiment 
and also can able to generate observable CP violation instead of small #13. Another thing we want 
to point out that negative value of (f> equally acceptable as far as it is small. It could not change 
mixing angles because their expressions depend on (f> through cos0 and sin 2 (j). Only Jcp and 5c p 
will change in sign which are unsettled according to the experiments. 

At the end we want to check whether the range of k and <ft can satisfy the double beta decay 
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Figure 7: Plot of Jcp with respect to k and cj) for the breaking of A4 in 13 element of in the 
unit of 10 -3 . We keep Am| and Am|i to their best fit values. 




Figure 8: Plot of 5qp with respect to k and for the breaking of A4 in "13" element of mo. We 
keep Am§ 2 and Am^i to their best fit values. 
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Figure 9: Plot of m e g = |(M iy ) ee | with respect to k and <p for the breaking of in "13" element 
of ijid ■ We keep Am^ and Am^ to their best fit values. 



bound m c fj = |(M„) ee | < 0.89 eV. In our model expression for this quantity is : 



i/wx , m o /9 + fc 2 (l + 2e) 2 + 6A;(l + 2e)cos0 . 

and it will be also only function of k and <f>. We plot this in Fig. [9] and it remains well below the 
experimental upper bound. 

Again we want to discuss about the mass pattern. Throughout our whole analysis in real as 
well as complex case we keep Arrig and Am^ tm to their best fit value and take the negative sign 
of Am 2 im . It corresponds to so called inverted ordering of neutrino mass. It is the feature near 
k = — 2. It is necessary to keep tuq > 0. Why we so fond of this region of k instead of region k = 1 
which can give the normal hierarchical mass spectrum. The reason is that the inverted ordering 
corresponds to the light neutrino mass scale mo — Am 2 im where mo — Am?, for normally ordered 
mass spectrum. So, from the point of view of observable CP violation, it is inevitable to choose 
larger value of mo because Jcp oc m(j- So, inverted hierarchical mass spectrum compatible with 
the observable CP violation. Now we extend our study of this model to leptogenesis. We want to 
see whether we can have appreciable leptogenesis compatible with baryon asymmetry for the same 
parameter space k and <fi after successful low energy data analysis for the feasible value of scale 
M . 
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4.2 Leptogenesis 



After successful predictions of low energy neutrino data we want to see whether this model can 
generate non-zero lepton-asymmetry with proper size and sign to describe baryon asymmetry. We 
keep the same right handed neutrino mass matrix Mr as before. The change only appear in Dirac 
type Yukawa coupling and hence in mo also. The diagonalisation of Mr gives 



M 



U J TB M R U^ B = diag{\A + De^\e i2a , \A\e rzA , \De l(p - A\e 



J2X 



(4.14) 



and hence the masses of the right handed neutrino are 



Mi = \A + De^\ 
M 3 

and the phases are 

tan 2a 



M \/l + k 2 + 2fecos. 



Mo 



\A\ = M 



\De l(p -A\= M Jl + k 2 - 2fccos 



k sin < 



1 + k cos 

The explicit form of diagonalising matrix is 

/ 

V = UtbUp = 



tan2A = tan 27 



ksm< 



\ 



k cos ( 







(4.15) 



(4.16) 



I / 





r 1 



AX 



(4.17) 







s 27 



where the expressions for the phases are given in Eq. (4.16). The mo matrix is no longer diagonal 



after explicit breaking of A4 symmetry. In the mass basis of right handed neutrino the modified 
Dirac mass term is m' D = m^V* . Hence the relevant matrix for describing Leptogenesis is : 



7 At ' 

h = m^m^) 



/ l-2e/3 
ee - i ( Q - A )/(3 v / 2) 



ee i(a-\) /(3V2) ee< a -"/)/y/3" 
l + 2e/3 



ee -i(A-7)/^6 



(4.18) 



To calculate lepton asymmetry as in Eq. (2.12) and Eq. ( 2.13| ) we need to calculate following 
quantities from matrix h: 



12) 



Im(/i 21 ) = —lm(h 

Im(/t| 2 ) = ~ lm ( h 23 
Ina(/i§ 1 ) = — Im(/i 



/V<sin2(a-A) /V< 



k sin ( 



72 



72 y/i + k 2 + 2kcos< 



/ 4 e 2 f^sin2(A — 7) / 4 e 2 ^ fcsin< 



13) 



12 
in 2 



fW uS in2(a- 7 ) _2fW u 



24 

_v 

12 



y/l + k 2 - 2k cos <f) 

A; sin 
< v /(l + A; 2 ) 2 -4fc 2 cos 2 l 



(4.19) 
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from Eq. (4.18) and taking Im(ftf) from Eq. (4.19) 



Calculating = Mf/Mf from Eq. 4.15, hi 
we calculate the self energy part of lepton asymmetry from Eq. (2.13) and vertex part of lepton 



asymmetry from Eq. (2.12 ). Adding both we obtain the following decay asymmetry of right handed 



neutrinos for all three generations 



El 



Mom^e ksincj) 
72irvl(l + k 2 + 2fecos ( 



1 + 2k 2 + 4fccos. 
k(k + 2 cos 4>) 



3(1 + k 2 + 6A:coS( 
A; cos (/> 



(2 + k 2 + 2/c cos 0) ln(2 + k 2 + 2k cos . 
1 + k 2 + 2k cos 6 



+ 



24(1 + A; 2 



1 + k 2 + 2k cos , 



x In 



2(1 + fc 2 



1 + k 2 - 2k cos < 



(4.20) 



£2 



Momoe 2 A;sin( 



x In 



727TO 2 



k 2 + 2k cos ( 



k(k + 2 cos ( 



1 3(1 

- + 



k + 2/ccos ( 



2 + k 2 + 2fecosi 
1 + k 2 + 2fccos( 



3(2 + k 1 



k(2 cos 
2k cos (/>) x In 



-k) 
2 + k 2 



(2 + k 2 + 2fccos( 
2fc cos 6 



1 + k 2 - 2fccos. 



(4.21) 



£3 



MoiriQe fcsini 



24vru 2 (l + A; 2 -2A;cos. 
2 + fc 2 - 2A;cosi 



6A; cos i 



1 + fc 2 - 2fcco Sl 



x ln(2 + fc' 



fecos « 
2k cos 0) — 



1 1 + 2A; 2 
- + 



4k cos ( 



k(k- 

8(1 + k 2 ) 
1 + k 2 -2k cos ( 



2 cos 
x In 



2(1 + /t 2 



1 + k 2 + 2A;cos( 



(4.22) 



CP asymmetry parameters £j are related to the leptonic asymmetry parameters through as 



til ~ n L 



E 



(4.23) 



where is the lepton number density, hl is the anti-lepton number density, s is the entropy 
density, m is the dilution factor for the CP asymmetry £j and g*i is the effective number of degrees 
of freedom [51] at temperature T = Mi. Value of g„ in the SM with three right handed Majorana 
neutrinos and one extra Higgs doublet is 116. The baryon asymmetry Yb produced through the 
sphaleron transmutation of Yl , while the quantum number B — L remains conserved, is given by 

EH 



Y B 



■Yl 



with 



8N F + 4N H 



(4.24) 



V" " 22N F + 13N H ' 

where Np is the number of fermion families and Njj is the number of Higgs doublets. The quantity 



8/23 in Eq. (4.24) for SM with two Higgs doublet. Now we introduce the relation between Yb 



18 



and i]b, where r\p is the baryon number density over photon number density n 7 . The relation is 

ESI 



Y B = 7.0394Yb, 



where the zero indicates present time. Now using the relations in Eqs. ( 4.23|4.24 4.25), 
and <7*i = 116 we have 



VB 



-3.23 x W- 2 J2^i 



This dilution factor approximately given by |50| I5"4"l |5"5] 

0.3 



with 



K, 



Hi 



(4.25) 
= 8/23 

(4.26) 
(4.27) 



where Tj is the decay width of iVj and Hi is Hubble constant at T = Mi. Their expressions are 



ha Mi 



and 



Hi = 



.M 2 



1 Attv 2 1 VW M P ' 

where v u = v sin/?, v = 246GeV and Mp = 1.22 x 10 19 GeV. Thus we have 

Mphu 



Ki 

For our model K\, K 2 and K3 are 

Mphu 



1.66 x ^Jg-vlMi 



(4.28) 



(4.29) 



K 2 
K 3 



(1 



26 ) 

3 > 



1.66 x 4vr^Tu2Mi 
M P h 22 



N x — x — 

X v 2 X (1 + k 2 +2/ccos0) 1 /2 



mn , 2e N 

, — = N x — jj- x (1 H 

1.66 x 4vr % /^w2M 2 ^ v 3 , 



Mp/l 33 



1 



1.66 x 4vr 



^ mo 

X «2 X (1 + fc2 _2jtcOS(/)) 1 /2 



(4.30) 



where TV 



1.66x4^^ 



5* = 9*1 ~ 5*2 ~ 5*3 ~ 116. Using Eq. (4.30) into Eq. (4.27) and from the 



expression of £j we can say that apart from logarithmic factor, K{ oc v^/tuq and £j oc mo/v 2 ,. So, 
baryon asymmetry will be independent of mo and i> u . They only appear through logarithmic factor 



in Ki- We consider tan/3 ~ 2.5. Substituting mo from Eq. (4.11), e from Eq. (4.10) and considering 



Mq some specific value into the expressions of £j, Kj we can have baryon asymmetry as function of k 



and 4> only. In Fig. 10 Fig. 11 and Fig. 12 we have plotted t]b as function of k and (j) in the unit of 
10~ 10 for three M values, M = 2 x 10 13 GeV, M = 5 x 10 13 GeV and M = 10 14 GeV rspectively. 
We have seen that the experimental value of rjp — 6.0 x 10~ 10 is obtainable in our model within 
the same range of k and <f> as in the low energy case. To see more explicitly the baryon asymmetry 
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m x 10 




Figure 10: Plot of baryon asymmetry t?b in unit of 10 with respect to k and <f> for the breaking 
of A4 in "13" element of uid ■ We keep Am| 2 and Am|i to their best fit values and have plotted 
for mass scale of right handed neutrino Mq = 2 x 10 13 GeV 



IB x 10 10 




Figure 11: Plot of baryon asymmetry 773 in unit of 10 with respect to k and <f> for the breaking 
of in "13" element of mp . We keep Am^ and Am^i to their best fit values and have plotted 
for mass scale of right handed neutrino Mq = 5 x 10 13 GeV 
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(/; x lO 10 




Figure 12: Plot of baryon asymmetry t/b in unit of 10 with respect to k and 4> for the breaking 
of A4 in "13" element of mo . We keep Am| 2 and Am^ to their best fit values and have plotted 
for mass scale of right handed neutrino Mq = 10 14 GeV 




Figure 13: We combine all three plots of baryon asymmetry for three right handed neutrino mass 
scale along with the WMAP value of baryon asymmetry 6 x 10~ 10 which is the plane surface. 
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Figure 14: Contour plot of baryon asymmetry tjb, 612 and, #23 in <^-fc plane for Am^ and Am| 
to their best fit values and for mass scale of right handed neutrino Mq = 2 x 10 13 GeV. 



plots we combine all three plots along with the observed baryon asymmetry value 6 x 10 10 which 



corresponds the plane surface in Fig. 13 The observed WMAP value of the baryon asymmetry 
curve intersect the lower curve (for Mq = 2 x 10 13 GeV) near the boundary of k and <j> variation. 
So, lower value of Mq could not generate observable baryon asymmetry. In the intersection region 
612 — 32° and 623 — 47°. If we allow that much of variation of #12 and #23 , we can have large low 
energy CP violation as well as baryon asymmetry with proper size and sign with Mq = 2 x 10 13 
GeV which is near the upper bound of right handed neutrino mass scale for generation of lepton as 
well as baryon asymmetry. If we more relax, we can easily see that the intersection of experimental 
and theoretical curve for Mq = 5 x 10 13 GeV and Mq = 10 14 GeV is in the lower value of k and cfr 
where well known neutrino mixing angles are more closer to their best fit value. 



Let us give a close look to the plot of t\b m Fig. 10 Fig. 11 and Fig. 12 near very low <fi and 
k ~ —2 region. This is the region where t]b become very large. From the expression of e\ and 
£2 it is clear that e\ and £2 both are singular at k + 2 cos <p = 0. This corresponds to equality of 
the masses M\ = M2 or x\2 = £21 = 1. This singularity can be avoided considering finite decay 
width of right handed neutrinos. We can able to maximize £\ and £2 and hence rjs using resonant 
condition Mi — M2 — Ti^/2. But, as we have already obtained the observed baryon asymmetry 
without resonance, it is not necessary to think about so finely tuned condition. Again in the region 
where the resonant condition is applicable, the Jcp is miserably small to observe through any 
experiments. 
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Figure 15: Contour plot of baryon asymmetry rj^, #12 and, 623 in cfr-k plane for Am^ and Am| 
to their best fit values and for mass scale of right handed neutrino Mq = 5 x 10 13 GeV. 
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Figure 16: Contour plot of baryon asymmetry 7713, O12 and, 623 in <p-k plane for Am^ and Am^ 
to their best fit values and for mass scale of right handed neutrino Mq = 10 14 GeV. 
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4.3 Final Analysis 



We end our analysis with the help of three contour plots of baryon asymmetry rjg = 6 x 10 10 
for three scales M = 2 x 10 13 GeV, M = 5 x 10 13 GeV and, M = 10 14 GeV in <j>-k plane. We 
insert the contours of 023 and 612 and manage to find the intersection of three contours for some 
reasonable value of #23 and #12. 

Case(I): Mq = 2 x 10 13 GeV, from Fig. 14 we are seeing that three contours 612 = 31.86°, 



#23 = 47° and, t]b = 6 x 10 10 are intersecting at a point ( 20.5, —2.082) in cp-k plane. So, the 
mixing angles are within nearly 2° variation about the best fit values. Obtained k, <f> value gives 
013 = 1.78°, J C p = 0.5 x 10~ 2 , 5 C p = 46°, m eff = 0.02423 eV, m = 0.05365 eV, e = -0.121 and 
fc + 2cos</> = -0.208. 



Case (II): Mq = 5 x 10 13 GeV, now from Fig. 15 we have the intersection of the contour 
Vb = 6 x 10~ 10 with the contours 6 12 = 33.1°, 23 = 46.2° at (16.0, -2.07) in cp-k plane. So, 6>i 2 
and 623 are more closer to their best fit values (nearly 1° deviation from their best fit value). At 
this point we have 0i 3 = 1°, J C p = 0.8 x 10~ 3 , 5 C p = 12°, m cff = 0.02175 eV, m = 0.0516 eV, 
e = -0.076 and k + 2 cos <h = -0.131. 



Case (III): M = 10 14 GeV, from the Fig. 16 the contours r] B = 6 x 10~ 10 , 9 l2 = 33.7° 



and, 023 = 46.0° have crossed in <p-k plane at (10.5, —2.08). So, higher scale of right handed 
neutrino mass helps to have very good value of $12 and 023- At the intersection point, we get 
013 = 0.8°, J C p = 0.23 x 10" 3 , 5 C p = 4.5°, m cS = 0.0194 eV, m = 0.051 eV, e = -0.065 and 
A; + 2 cos = -0.11. 

The small value of e compatible with all experimental results. So, we can demand that A4 is 
an approximate symmetry. 



4.4 Corelation of CP violating phases in this model 

In this model everything is determinable in terms of parameter k and (j). Well measured quantities 
fix the value of those parameters. So, value of the rest of the physical quantities (some of them 
not so well measured in experiment like 0i3, m e g and some of them yet to measure in experiment 
like Jcp, 5q p and the Majorana phases also) are obtainable in this model. Question may arise 
whether we can have any relations among the phases in this model or not. First of all, there are two 
kinds phases, low energy and high energy phases. High energy phases are responsible to generate 
the lepton asymmetry. The low energy phases are responsible for determining low energy leptonic 
CP violation. Low energy phases are in two type, one is the lepton no preserving CP violating 
phase 5c p and another two are the lepton no breaking CP violating phases. In general, all phases 
are independent, meaning that there are no correlation among the phases between high and low 
energy sector, and also phases inside a particular sector are not correlated. For three generations 
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of neutrinos, there are three key phases which are responsible for leptogenesis. Those are phases in 
h\ 2 , ^13) and h^. From matrix h given in Eq. (4.18), we have 



$1 = arg(>f 2 ) = 2a 
It leads to the relation, 



$ 2 = arg(/i? 3 ) = 2(q - 7) 



$2 = $1 +$3- 



<!>■: 



arg(/i 



23 ) 



-2 7 . 



(4.31) 



(4.32) 



For a given k and the a and 7 are known from Eq. (4.16). Hence <J?i, $2 and $3 are individually 
determinable phases. But, in this model values of those high energy phases follow the relation given 
in Eq. Q4.32I). 



Now, let us give a fresh look to the leptonic mixing matrix. In Eq. (4.3) we have given the 
tri-bimaximal rotated form of the neutrino mass matrix. Keeping terms upto first order in e we 
obtain the diagonalising matrix in the following form, 



V = U TB 



( 



1 



-eXe- i0x 
eXe i0x 1 
eYe idY eXe i6z 



-eYe~ l0Y 
-eZe~ i0 z 
1 



\ 



/ 



(4.33) 



where the X, Y, Z and the associated phases are completely known function of k and (p. An 



additional phase matrix Vp is needed to make masses of light neutrino real, from Eq. (4.3) we have 
the phase matrix 



V P 



( e i{n/2-a) o o 

e™l 2 
\ e ^/ 2 -^ 



\ 



(4.34) 



To obtain the CKM form of mixing matrix we need to rotate V by two diagonal phase matrix, let 
U P1 = diag(e iQl , e ic * 2 , e iQ3 ) and U P2 = diag(e i/?1 , e ih , e^ 3 ). So, we have 



V = U Pl {U P1 V'U P2 } U P2 = U p1 VckmU P 2 

Now the with #13 small we can write 

{V C km)ii = co S 6 12 = Vl^-M = IV^l^^-M 
(V C km)i2 = sinfl 12 = V{ 2 e< a ^M = \V[ 2 \e l ^ +a ^ 



(4.35) 



{Vckm)v& = sm0i 3 e 



-iSci 



V[ e^ ai ~^ = \V' \ e i ^ Sl+ai ~^ 3 \ 



(4.36) 



and more six relations. But these three are sufficient for our discussions. The phases associated 
to V' elements, like ipi, 4>i, and 5\ associated to V[i, V{ 2 and V{ 3 respectively, are completely 
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determinable in terms of k, 4> using functional form of e, X, Y, Z and their associates phases. Now 
from the Eq. (4.36) we obtain the following phase relations 



tpi + ax - fix = 0, 4>\ + ai - /3 2 = 0, 5i + a\ - /3 3 = -S C p. 
Now the form of total mixing matrix is, 

/ e i(7r/2-a+/3i) 



(4.37) 



V 



V'Vp = U ] p1 V C kmUp 2 Vp = U ] p1 V C km 

( e i(A-/3 3 +7-") 



V 



o 
o 



{e^/ 2 -^ + ^U P1 }Vc 



CRM 









e i(7r/2+/3 2 ) 








= i(7r/2- 7 +/3 3 ) 



V 



o 





e i(7+/32-/3 3 ) o 

1 



(4.38) 



This phase part in the parenthesis can be absorbed to charged lepton fields and the remaining part 
gives the leptonic mixing matrix of the form Vpmns = Vckm x diag(e* QM , e*^ A/ , 1), where the 
cum and (5m are the two Majorana phases of leptonic mixing matrix. From Eq. (4.38) and using 
relations in Eq. (4.37) and Eq. (4.31), we have the Majorana phases, 

$2 



«M 



-5cp 



M 



Scp~y + ® 



(4.39) 



where ^ = ipi — 5± and $ = 4>\ — S± are known function of k and (f>. In Eq. (4.39) we have the 
relation of low energy and high energy phases. So, in our model we have correlation among the CP 
violating phases. 



5 Conclusion 



We have shown that non-zero U e ^ is generated in a softly broken symmetric model through see- 
saw mechanism incorporating single parameter perturbation in mp> in single element. First, we have 
studied all possible nine cases to explore the mixing angles considering all model parameters real. 
The extent of 9\s investigated, keeping the experimental values of present solar and atmospheric 
mixing angles. Among all nine possible texture of mp some of them generates non-zero 613. Out 
of those non-zero #13 generating textures of mo we find that breaking at '12' and '13' elements are 
encompassing the best values of Q\i and #23- However, the reach of #13 in those cases are around 
1°. Considering one of the parameter complex we extend our analysis with one of the most suitable 
texture of mp> with breaking at '13' element. We have calculated mixing angles and neutrino mass 
squared differences in terms four model parameters (mo, e, k, 4>). We restrict model parameters 
utilising the well measured quantities Airig, Am^ (m , #12 and 6*23 and we have obtained #13 (upto 2°) 
and large Jcp (— 10 _2 )and |(M„) ee | well below the present experimental upper bound. In addition 
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to that a large Sep is also obtained. Further study on leptogenesis is also done and the present 
WMAP value of baryon asymmetry is obtained for a right handed neutrino mass scale Mq ~ 10 13 
GeV. In our model, we have seen that the phases responsible for the leptogenesis are correlated. We 
also find out the relations among low energy CP violating phases and the lepton asymmetry phases. 
Small A4 symmetry breaking parameter e, is sufficient to describe the all low energy neutrino data 
and high energy CP violation (leptogenesis). So, A4 symmetry is an approximate symmetry. 



Here we consider breaking of A4 symmetry in all other entries of mo- Case (i) is already discussed 
in the text. 

(ii) Breaking at '11' element : In this case mo is given by 

fv u 



m D 




The mass eigenvalues are 
mi = -2jDTA) 



1 + 



4e 



2A 



m 2 = — 



1 + 



2e 



m 3 



2(D — A) 



and the three mixing angles come out as 
. 1 2(2 + k) 

sm0i2= 7i + ^r 



sin tf 23 = - 



V2 



sin = 



The solar and atmospheric mass differences and their ratio are 



Am ° = WTW ^ Kk + 2)+ Hk2 + 2k ~ 1] 



Am 2 - m ° 
atm ~ 3(1 - k) 2 



Am; 



[3fc(2 — k) - 4e(fc - 1) 



(A.l) 



(A.2) 



(A.3) 



- (k - l) 2 [3k(k + 2) + 4e(fc 2 + 2k - 1)] 

" AmL " (k + l) 2 [3k(2 -k)- 4e(k - l) 2 ] 

The obtained expression for e in terms of model parameter k and experimentally known R: 

3k [R{2-k){k + l) 2 -jk-l) 2 {k + 2)] 
€ ~ 4(k-l) 2 [R(k + l) 2 + k 2 + 2k - 1] 



(A.4) 



(A.5) 
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iii) Breaking at '33' element : In this case, the structure of tod is given by 

(\ \ 



m D 



fVu 



1 
\0 1 + eJ 



(A.6) 



Mass eigenvalues are 



mi = — 



2(D + A) 



1 + 



m 3 = - 



2{D - A) 



+ 



and the angels are 

. . 1 2+fc 



sin C/i2 



V3 3y/3k 



e sin y 2 3 = - 



1 fc(4 - fc) 

V2 £ 6 v / 2(2-/c) 



sin 



13 



efc(l - fe) 
3^(2 - fc) 



(A.7) 



(A.8) 



The solar and atmospheric mass differences and their ratio are 



Am 



rrin 







Am 



3(1 + A;) 2 



atm 



3(1 -fc) 2 



3/c(fc + 2) + 2e(2k 2 + Ak + 1) 



3k(2 - fc) - 2e(2/c 2 - 4k - 1) 



R 



Am 2 Q _( k ~ l ) 2 [ 3k ( k + 2) + 2e(2A; 2 + Ak + 1)] 



^ m atm 



(k + l) 2 [3fc(2 - k) - 2e(2k 2 - Ak - 1)] 

The obtained expression for e in terms of model parameter k and experimentally known R: 

_ 3k [R(2- k){k + l) 2 - jk-l) 2 {k + 2)] 

6 ~ ~2 p - l) 2 (2fc 2 + Ak + 1) + i?(£; + 1) 2 (2A: 2 - 4fc - 1)] 

iv) Breaking at '12' element : In this case, the structure of mo is given by 

(\ e 0\ 



(A.9) 



(A.10) 



mo 



fVu 



Mass eigenvalues are 



mi 



2{D + A) 



2e 
~3 



m 2 



1 

Vo o 1/ 



2A \ ^ 3 



(A.ll) 



m 3 



2(Z>- A) 



(A.12) 



and the three mixing angles come out as 



sin &i2 



1 2k + 1 



V3 SV3k 



e sin t>23 = - 



k(l - k) 
^/2 "6\/2(2 — k) 



— € 



sin 013 = - 



e k 2 -k-3 
3^2 (k - 2) 



(A.13) 
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The solar and atmospheric mass differences and their ratio are 



m 



Am 



3(1 + k) 2 



atm 



3k(k + 2) + 4e(/c 2 + 2k + 2) 



3fc(2 - fc) -4e(l - A;) 5 



3(1 - A;) 2 

(A; - l) 2 [3ife(fc + 2) + 4e(/c 2 + 2A; + 2)] 



(fc + l) 2 [3k(2 - k) - 4e(l - A;) 2 ] 

The obtained expression for e in terms of model parameter k and experimentally known R: 

_ 3k [R{2-k){k + l) 2 -jk-l) 2 {k + 2)] 
e ~ 4(1- k) 2 [(At 2 + 2k + 2) + R(k + l) 2 ] 

v) Breaking at '13' element : In this case, the structure of mp is given by 



(A.14) 



(A.15) 



m D 




1 





e\ 





1 











J 



The mass eigenvalues are 

f 2 vl (, 2 £ 
mi = -2(DTA){ 1 -J 



m * = -2t 1 + y 



m 3 = - 



2(D- A) 



and the three mixing angles come out as 

1 2A: + 1 

sin 6»i3 



sm #23 = - 
e A; 2 -A; -3 



1 +e fc(l-fc) 



6^(2 - k) 



3^/2 (k — 2) 

The solar and atmospheric mass differences and their ratio are 

2 

Am| = m ° [3k(k + 2) + 4e(A; 2 + 2A; + 2) 



Ara <,(m = 



[3fe(2 - fc) - 4e(l - A;) 5 



3(1 - A;) 2 

(A; - l) 2 [3k(k + 2) + 4e(A: 2 + 2A; + 2)] 



(A.16) 



(A.17) 



(A.18) 



(A; + l) 2 [3A?(2 - k) - 4e(l - A:) 2 ] 
The obtained expression for e in terms of model parameter k and experimentally known R 



3k [R(2-k)(k + l) 2 - (k-l) 2 (k + 2)] 
~ 4(1 - A:) 2 [(k 2 + 2 A; + 2) + R(k + l) 2 ] 



(A.19) 



(A.20) 
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vi) Breaking at '23' element : In this case, the structure of mo is given by 

/l 0\ 
m D = — 1 6 
\0 1/ 



(A.21) 



The mass eigenvalues are 
mi ~ 



■III) - A) V.' ' 3 



m 2 



V 3 



m 3 



2(L> - A) 



(1-6) 



and the three mixing angles come out as 
1 e 2 + fc 

smfc'12 



3y/3 fc 



sm t/ 2 3 



1 



V2 2^2 



sin #13 = 



(A.22) 



(A.23) 



The solar and atmospheric mass differences and their ratio are 



Ami 



rrin 



3(1 + A;) 2 
,2 



[3fc(fc + 2) + 2e(2fc 2 + 4fc + 1) 



Am ' tm = 3(1-^)2 [ 3fc ( 2 " fc ) " 2e ( 2fc2 " 4fc + 5 ) 

Am| _ (fc - l) 2 [3fc(fc + 2) + 2e(2fc 2 + 4fc + 1)] 
~ Am 2 im ~ (fc + l) 2 [3fc(2 - k) - 2e(2k 2 - 4k + 5)] 



(A.24) 



The obtained expression for e in terms of model parameter k and experimentally known R: 

_ 3k [R(2- k){k + l) 2 - {k- l) 2 (fc + 2)] 

6 ~ Y [R(l + A;) 2 (2fc 2 - 4fc + 5) + (k - l) 2 (2fc 2 + 4k + 1)] 

vii) Breaking at '21' element : In this case, the structure of mc is given by 

/1 (r 

mD = ~Ji e 1 
Vo 1 



(A.25) 



(A.26) 



The mass eigenvalues are 
mi = -2jDTA) 



2e 



m 2 = — 



2A V 3 



m 3 



2(D — A) 



(A.27) 



and the three mixing angles come out as 

1 e 1-k 

sin 



sin 6*23 = 
e (3 + fc 2 



1 e fc(fc-l) 
72 + QV2 (2 - fc) 



4k) 



3V2 (2 - fc) 



(A.28) 
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The solar and atmospheric mass differences and their ratio are 

2 

mf 



Am 2 G = 3 ( 1 +° fc )2 [ 3k( - k + 2 ) + 4e ( fc2 + 2k + 2 ) 



Am 



atm 



3(1 - 



[3ife(2 - k) 



4e(l - A;) 2 



Am| (A:-l) 2 [3fc(£; + 2) + 4e(fc 2 + 2fc + 2)] 



(fc + 1) 2 [3fc(2 - fc) - 4e(l - fc) 2 ] 

The obtained expression for e in terms of model parameter A; and experimentally known i?: 

3fc [i?(2-fc)(fc + l) 2 -(fc + 2)(fc-l) 2 ] 
6 ~ 4(1 - /c) 2 [R(l + k) 2 + k 2 + 2k + 2] 

viii) Breaking at '31' element : In this case, the structure of mr> is given by 



m D 




The mass eigenvalues are 
mi = -2(DTA) 



2e 
3~ 



2„,2 



1712 



2A 



1 + 



2c 



W3 



2{D — A) 



and the three mixing angles come out as 

1 e 1-k 

sin 0i3 



sin #23 



e 3 + Ar 



1 e fc(l-ife) 

7^ + qV^ (2 - fc) 



4A; 



3^ 2-fe 

The solar and atmospheric mass differences and their ratio are 



a 2 m n 

Am « = 3(iT^ 



Am 



m- r 



atm 



3(1 - k) 2 



3k(k + 2) + 4e(/c 2 + 2fc + 2) 



3fc(2 - fc) -4e(l - kf 



Am| (fc-l) 2 [3fc(/c + 2) + 4e(/c 2 + 2A; + 2)] 



Am ltm 



(A.29) 



(A.30) 



(A.31) 



(A.32) 



(A.33) 



(A; + l) 2 [3/c(2 — fc) — 4e(l - k) 2 } 

The obtained expression for e in terms of model parameter k and experimentally known R 

3k [R{2-k){k + l) 2 -{k-l) 2 {k + 2)] 
6 ~ 4(1 - k) 2 [R(l + k) 2 + k 2 + 2k + 2] 



(A.34) 



(A.35) 
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ix) Breaking at '32' element : In this case, the structure of mo is given by 

1 (T 

fv u I 

m D 




The mass eigenvalues are 



mi 



2(D + A) V 1+ 3 



m 2 



2A I 3 



m 3 



f 2 i> 2 



2(L> - A) 



(1-6) 



and the three mixing angles come out as 

1 e 2 + k 

/3 



sm 12 



sin 6/ 2 3 



1 e 

/2 2^2 



3V3 

The solar and atmospheric mass differences and their ratio are 



sin #13 = 



Ami 



Ami,, 



mf, 



3(1 + k) 2 



mZ 



atm 



3k(k + 2) + 2e(2k 2 + 4k + 1] 



3k(2 - k) - 2e(2k 2 - 4k + 5) 



R 



Ami 



3(1 - k) 2 

(k - l) 2 [3k{k + 2) + 2e{2k 2 + 4k + 1)] 



(A.36) 



(A.37) 



(A.38) 



(k + l) 2 [3fc(2 - k) - 2e(2k 2 - 4k + 5)] 

The obtained expression for e in terms of model parameter k and experimentally known R: 

_ 3k [R{2-k){k + l) 2 - jk-l) 2 {k + 2)] 

6 ~ ~2 [R(l + k) 2 (2k 2 - 4k + 5) + (1 - k) 2 (2k 2 + 4k + 1)] 



(A.39) 



(A.40) 
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